SLAC-PUB-13641 
BAB4R-PUB-09/006 



Measurement of the 77* — > n transition form factor 

B. Aubcrt, Y. Karyotakis, J. P. Lees, V. Poireau, E. Prencipe, X. Prudent, and V. Tisserand 
Laboratoire d'Annecy-le-Vieux de Physique des Particules (LAPP), 
Universite de Savoie, CNRS/IN2P3, F-74941 Annecy-Le-Vieux, France 

J. Garra Tico and E. Granges 

Umversitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain 

M. Martinelli°^ A. Palano"^ and M. Pappagallo"'' 
INFN Sezione di Barf^ ; Dipartimento di Fisica, Universitd di Barf , 1-70126 Bari, Italy 

G. Eigen, B. Stugu, and L. Sun 
University of Bergen, Institute of Physics, N-5007 Bergen, Norway 

M. Battaglia, D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, L L. Osipenkov, K. Tackmann, and T. Tanabe 
Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 

C. M. Hawkes, N. Soni, and A. T. Watson 
University of Birmingham, Birmingham, B15 2TT, United Kingdom 

H. Koch and T. Schroeder 
Ruhr Universitdt Bochum, Institut fiir Experimentalphysik 1, D-44780 Bochum, Germany 

D. J. Asgcirsson, B. G. Fulsom, C. Hearty, T. S. Mattison, and J. A. McKenna 
University of British Columbia, Vancouver, British Columbia, Canada V6T IZl 

M. Barrett, A. Khan, and A. Randle-Conde 

Brunei University, Uxbridge, Middlesex UBS 3PH, United Kingdom 

V. E. Bhnov, A. D. Bukin,* A. R. Buzykaev, V. P. Druzhinin, V. B. Gohibev, 
A. P. Onuchin, S. L Serednyakov, Yu. L Skovpen, E. P. Solodov, and K. Yu. Todyshev 
Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

M. Bondioh, S. Curry, L Eschrich, D. Kirkby, A. J. Lankford, P. Lund, M. Mandelkern, E. C. Martin, and D. P. Stoker 

University of California at Irvine, Irvine, California 92697, USA 

H. Atmacan, J. W. Gary, F. Liu, O. Long, G. M. Vitug, Z. Yasin, and L. Zhang 
University of California at Riverside, Riverside, California 92521, USA 

V. Sharma 

University of California at San Diego, La Jolla, California 92093, USA 

C. Campagnari, T. M. Hong, D. Kovalskyi, M. A. Mazur, and J. D. Richman 

University of California at Santa Barbara, Santa Barbara, California 93106, USA 

T. W. Beck, A. M. Eisner, C. A. Heusch, J. Kroseberg, W. S. Lockman, 
A. J. Martinez, T. Schalk, B. A. Schumm, A. Seiden, L. Wang, and L. O. Winstrom 
University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 

C. H. Cheng, D. A. DoU, B. Echenard, F. Fang, D. G. Hithn, L Narsky, T. Piatenko, and F. C. Porter 
California Institute of Technology, Pasadena, California 91125, USA 

R. Andreassen, G. MancineUi, B. T. Meadows, K. Mishra, and M. D. Sokoloff 



2 



University of Cincinnati, Cincinnati, Ohio 45221, USA 

P. C. Bloom, W. T. Ford, A. Gaz, J. F. Hirschauer, M. Nagel, U. Nauenberg, J. G. Smith, and S. R. Wagner 

University of Colorado, Boulder, Colorado 80309, USA 

R. Ayad,'f W. H. Toki, and R. J. Wilson 

Colorado State University, Fort Collins, Colorado 80523, USA 

E. Feltresi, A. Hauke, H. Jasper, T. M. Karbach, J. Merkel, A. Petzold, B. Spaan, and K. Wacker 
Technische Universitdt Dortmund, Fakultdt Physik, D-44221 Dortmund, Germany 

M. J. Kobel, R. Nogowski, K. R. Schubert, R. Schwierz, and A. Volk 
Technische Universitdt Dresden, Institut fiir Kern- und Teilchenphysik, D-01062 Dresden, Germany 

D. Bernard, E. Latour, and M. Vcrdcri 
Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Polytechnique, F-91128 Palaiseau, France 

P. J. Clark, S. Playfer, and J. E. Watson 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

M. A^dreotti'^^ D. Bettoni", C. Bozzi'', R. Calabrese''^ A. Cecchi''^ G. Cibinetto"^ E. Fioravanti'^^ 

P. Franchi^i''^ E. Luppi"^ M. Munerato°^ M. Negrini°^ A. Petrella"^ L. Piemontese", and V. Santoro'''' 
INFN Sezione di Ferrara"^ ; Dipartimento di Fisica, Universitd di Ferrara^ , 1-44100 Ferrara, Italy 

R. Baldini-Fcrroli, A. Calcaterra, R. de Sangro, G. Finocchiaro, 
S. Pacetti, P. Patteri, I. M. Peruzzi,* M. Piccolo, M. Rama, and A. Zallo 
INFN Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 

R. Contri°^ E. Guido, M. Lo Vetere"^ M. R. Monge'^^ S. Passaggio", C. Patrignani''^ E. Robutti", and S. Tosi"^ 
INFN Sezione di Genova"' ; Dipartimento di Fisica, Universitd di Genova'' , I-I6I46 Genova, Italy 

K. S. Chaisanguanthum and M. Morii 

Harvard University, Cambridge, Massachusetts 02138, USA 

A. Adametz, J. Marks, S. Schenk, and U. Uwer 
Universitdt Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 

F. U. Bernlochner, V. Klose, and H. M. Lacker 
Humboldt-Universitdt zu Berlin, Institut fiir Physik, Newtonstr. 15, D-12489 Berlin, Germany 

D. J. Bard, P. D. Daunccy, and M. Tibbetts 
Imperial College London, London, SW7 2AZ, United Kingdom 

P. K. Behera, M. J. Charles, and U. Mallik 

University of Iowa, Iowa City, Iowa 52242, USA 

J. Cochran, H. B. Crawley, L. Dong, V. Eyges, W. T. Meyer, S. Prell, E. I. Rosenberg, and A. E. Rubin 

Iowa State University, Ames, Iowa 5001 1-3160, USA 

Y. Y. Gao, A. V. Gritsan, and Z. J. Guo 
Johns Hopkins University, Baltimore, Maryland 21218, USA 

N. Arnaud, J. Bequilleux, A. D'Orazio, M. Davier, D. Derkach, J. Firmino da Costa, 
G. Grosdidier, F. Le Dibcrdcr, V. Lepeltier, A. M. Lutz, B. Malaescu, S. Pruvot, 
P. Roudeau, M. H. Schune, J. Serrano, V. Sordini, § A. Stocchi, and G. Wormser 
Laboratoire de I'Accelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, 
Centre Scientifique d'Orsay, B. P. 34, F-91898 Orsay Cedex, France 



3 



D. J. Lange and D. M. Wright 
Lawrence Livermore National Laboratory, Livermore, California 94550, USA 

I. Bingham, J. P. Burke, C. A. Chavez, J. R. Fry, E. Gabathuler, 
R. Garnet, D. E. Hutchcroft, D. J. Payne, and C. Touramanis 

University of Liverpool, Liverpool L69 7ZE, United Kingdom 

A. J. Bevan, C. K. Clarke, F. Di Lodovico, R. Sacco, and M. Sigamani 

Queen Mary, University of London, London, El 4NS, United Kingdom 

G. Cowan, S. Paramcsvaran, and A. C. Wren 
University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 

D. N. Brown and C. L. Davis 
University of Louisville, Louisville, Kentucky 40292, USA 

A. G. Dcnig, M. Fritsch, W. Gradl, and A. Hafner 
Johannes Cutenberg-Universitdt Mainz, Institut fiir Kernphysik, D-55099 Mainz, Cermany 

K. E. Alwyn, D. Bailey, R. J. Barlow, G. Jackson, G. D. Lafferty, T. J. West, and J. I. Yi 
University of Manchester, Manchester M13 9PL, United Kingdom 

J. Anderson, C. Chen, A. Jawahery, D. A. Roberts, G. Simi, and J. M. Tuggle 
University of Maryland, College Park, Maryland 20742, USA 

C. Dallapiccola, E. Salvati, and S. Saremi 

University of Massachusetts, Amherst, Massachusetts 01003, USA 

R. Cowan, D. Dujmic, P. H. Fisher, S. W. Henderson, G. ScioUa, M. Spitznagel, R. K. Yamamoto, and M. Zhao 
Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 

P. M. Patel, S. H. Robertson, and M. Schram 
McGill University, Montreal, Quebec, Canada H3A 2T8 

A. Lazzaro"^ V. Lombardo", F. Palombo''^ and S. Stracka"'' 
INFN Sezione di Milano"' ; Dipartimento di Fisica, Universitd di Milano^ , L20133 Milano, Italy 

J. M. Bauer, L. Cremaldi, R. Godang,^ R. Kroeger, P. Sonnek, D. J. Summers, and H. W. Zhao 
University of Mississippi, University, Mississippi 38677, USA 

M. Simard and P. Taras 
Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 3J7 

H. Nicholson 

Mount Holyoke College, South Hadley, Massachusetts 01075, USA 

G. De Nardo"^ L. Lista", D. Monorchio"^ G. Onorato''^ and C. Sciacca"^ 
INFN Sezione di NapoW^ ; Dipartimento di Scienze Fisiche, 
Universitd di Napoli Federico if, 1-80126 Napoli, Italy 

G. Raven and H. L. Snoek 
NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

C. P. Jessop, K. J. Knoepfel, J. M. LoSecco, and W. F. Wang 
University of Notre Dame, Notre Dame, Indiana 46556, USA 

L. A. Corwin, K. Honscheid, H. Kagan, R. Kass, J. P. Morris, 
A. M. Rahimi, J. J. Regensburger, S. J. Sekula, and Q. K. Wong 



4 



Ohio State University, Columbus, Ohio 43210, USA 

N. L. Blount, J. Brau, R. Frey, O. Igonkina, J. A. Kolb, M. Lu, 
R. Rahmat. N. B. Sincv, D. Strom, J. Strubc, and E. Torrence 
University of Oregon, Eugene, Oregon 97403, USA 

G. Castelli''\ N. Gagliardi"^ M. Margoni"^ M. Morandin'^, 
M. Posocco", M. Rotondo'', F. Simonetto''^ R. Stroili"^ and C. Voci'^'' 
INFN Sezione di Padova^ ; Dipartimento di Fisica, Universitd di Padovci' , 1-35131 Padova, Italy 

P. del Amo Sanchez, E. Ben-Haim, G. R. Bonneaud, H. Briand, J. Chauveau, 
O. Hamon, Ph. Leruste, G. Marchiori, J. Ocariz, A. Perez, J. Prendki, and S. Sitt 
Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, Universite Pierre et Marie Curie-Paris6, 
Umversite Dems Diderot- Paris!, F-75252 Paris, France 

L. Gladney 

University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 

M. Biasini"'' and E. Manoni"^ 
INFN Sezione di Perugia"^ ; Dipartimento di Fisica, Universitd di Perugia'' , 1-06100 Perugia, Italy 

C. Angelini''^ G. Batignani"^ S. Bettarini'^^ G. Calderini'^^** M. Carpinelli"^tt A. Cervelli'^^ F. Forti"^ 
M. A. Giorgi''^ A. Lusiani'''=, M. Morganti°^ N. Neri"^ E. Paoloni''^ G. Rizzo"^ and J. J. Walsh'^ 
INFN Sezione di Pisa"^ ; Dipartimento di Fisica, Universitd di Pisa'' ; Scuola Normale Superiore di Pisa'', 1-56127 Pisa, Italy 

D. Lopes Pegna, C. Lu, J. Olsen, A. J. S. Smith, and A. V. Telnov 
Princeton University, Princeton, New Jersey 08544, USA 

F. AnuUi'^, E. Baracchini''^ G. Cavoto", R. Faccini'^^ F. Ferrarotto", F. Ferroni'^^ M. Gaspero'^^ 
P. D. Jackson", L. Li Gioi°, M. A. Mazzoni'^, S. Morganti'^, G. Piredda", F. Renga"^ and C. Voena" 

INFN Sezione di Roma"; Dipartimento di Fisica, 
Universitd di Roma La Sapienzc^ , 1-00185 Roma, Italy 

M. Ebert, T. Hartmann, H. Schroder, and R. Waldi 
Universitdt Rostock, D-18051 Rostock, Germany 

T. Adye, B. Franek, E. O. Olaiya, and F. F. Wilson 

Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

S. Emery, L. Esteve, G. Hamel de Monchenault, W. Kozanecki, G. Vasseur, Ch. Ycche, and M. Zito 
CEA, Irfu, SPP, Centre de Saclay, F-91191 Gif-sur-Yvette, France 

M. T. Allen, D. Aston, R. Bartoldus, J. F. Benitez, R. Cenci, J. P. Coleman, M. R. Convery, 
J. C. Dingfelder, J. Dorfan, G. P. Dubois-Felsmann, W. Dunwoodie, R. C. Field, M. Franco Sevilla, 
A. M. Gabareen, M. T. Graham, P. Grenier, C. Hast, W. R. Innes, J. Kaminski, M. H. Kelsey, H. Kim, 

P. Kim, M. L. Kocian, D. W. G. S. Leith, S. Li, B. Lindquist, S. Luitz, V. Luth, H. L. Lynch, 
D. B. MacFarlane, H. Marsiske, R. Messner,* D. R. MuUer, H. Neal, S. Nelson, C. P. O'Grady, L Oftc, 
M. Perl, B. N. RatclifF, A. Roodman, A. A. Salnikov, R. H. Schindler, J. Schwiening, A. Snyder, D. Su, 
M. K. Sullivan, K. Suzuki, S. K. Swain, J. M. Thompson, J. Va'vra, A. P. Wagner, M. Weaver, C. A. West, 
W. J. Wisniewski, M. Wittgen, D. H. Wright, H. W. Wulsin, A. K. Yarritu, C. C. Young, and V. Ziegler 
SLAC National Accelerator Laboratory, Stanford, California 94309 USA 

X. R. Chen, H. Liu, W. Park, M. V. Purohit, R. M. White, and J. R. Wilson 
University of South Carolina, Columbia, South Carolina 29208, USA 



P. R. Burchat, A. J. Edwards, and T. S. Miyashita 



5 

Stanford University, Stanford, California 94305-4060, USA 

S. Ahmed, M. S. Alam, J. A. Ernst, B. Pan, M. A. Saeed, and S. B. Zain 

State University of New York, Albany, New York 12222, USA 

A. Soffer 

Tel Aviv University, School of Physics and Astronomy, Tel Aviv, 69978, Israel 

S. M. Spanier and B. J. Wogsland 
University of Tennessee, Knoxville, Tennessee 37996, USA 

R. Eckmann, J. L. Ritchie, A. M. Ruland, C. J. Schilling, R. F. Schwitters, and B. C. Wray 
University of Texas at Austin, Austin, Texas 78712, USA 

B. W. Dmmmond, J. M. Izcn, and X. C. Lou 
University of Texas at Dallas, Richardson, Texas 75083, USA 

F. Bianchi''^ D. Gamba°^ and M. PcUiccioni'''' 
INFN Sezione di Torino"" ; Dipartimento di Fisica Sperimentale, Universitd di Torino'' , 1-10125 Torino, Italy 

M. Bomben"^ L. Bosisio°^ C. Cartaro"\ G. Delia Ricca''^ L. Lanceri'^^ and L. Vitale"'' 
INFN Sezione di Trieste"" ; Dipartimento di Fisica, Universitd di Trieste*" , 1-34127 Trieste, Italy 

V. Azzolini, N. Lopez-March, F. Martinez- Vidal, D. A. Milanes, and A. Oyanguren 
IFIC, Universitat de Valencia- C SIC, E-46071 Valencia, Spain 

J. Albert, Sw. Banerjee, B. Bhuyan, H. H. F. Choi, K. Hamano, G. J. King, 
R. Kowalewski, M. J. Lewczuk, L M. Nugent, J. M. Roney, and R. J. Sobie 

University of Victoria, Victoria, British Columbia, Canada V8W 3P6 

T. J. Gershon, P. F. Harrison, J. Ihc, T. E. Latham, G. B. Mohanty, and E. M. T. Puccio 
Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 

H. R. Band, X. Chen, S. Dasu, K. T. Flood, Y. Pan, R. Prepost, C. O. Vuosalo, and S. L. Wu 
University of Wisconsin, Madison, Wisconsin 53706, USA 

We study the reaction e^e^ e^e^vr" and measure the 77* tt" transition form factor in 
the momentum transfer range from 4 to 40 GeV'^. The analysis is based on 442 fb^^ of integrated 
luminosity collected at PEP-II with the BABAR detector at e'^e~ center-of-mass energies near 10.6 
GeV. 

PACS numbers: 14.40.Aq, 13.40.Gp, 12.38.Qk 



I. INTRODUCTION 



'Deceased 

tNow at Temple University, Philadelphia, Pennsylvania 19122, 
USA 

■tAlso with Universita di Perugia, Dipartimento di Fisica, Perugia, 
Italy 

§Also with Universita di Roma La Sapienza, 1-00185 Roma, Italy 
^Now at University of South Alabama, Mobile, Alabama 36688, 
USA 

**Also with Laboratoire de Physique Nucleaire et de Hautes Ener- 
gies, IN2P3/CNRS, Universite Pierre et Marie Curie-Paris6, Uni- 
versite Denis Diderot-Paris7, F-75252 Paris, France 
ttAlso with Universita di Sassari, Sassari, Italy 



In this paper we study the process 

e+e- ^ e+e-7r°, (1) 

where the final state 7r° is produced via the two-photon 
production mechanism illustrated by Fig. 1. We measure 
the differential cross section for this process in the single 
tag mode where one of the outgoing electrons^ (tagged) is 
detected while the other electron (untagged) is scattered 



^ Unless otherwise specified, we use the term "electron" for either 
an electron or a positron. 
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FIG. 1: The Feynman diagram for the e"^e e^e tv'^ two- 
photon production process. 



at a small angle. The ir^ is observed through its decay 
into two photons. The tagged electron emits a highly off- 
shell photon with the momentum transfer q\ = —Q^ = 
{p — p'Y' , where p and p' are the four momenta of the 
initial and final electrons. The momentum transfer to 
the untagged electron is near zero. The differential cross 
section for pseudoscalar meson production da{e^e'' 
e^e~'K^)/dQ^ depends on only one form factor, F{Q^), 
which describes the 77* — > 7r° transition. To relate the 
differential cross section to the transition form factor we 
use the formulae for the e+e^ — > e^e^Tr" cross section in 
Eqs. (2.1) and (4.5) of Rcf. [1]. 

At large momentum transfer, Q^, perturbative QCD 
(pQCD) predicts that the transition form factor can be 
represented as a convolution of a calculable hard scatter- 
ing amplitude for 77* qq with a nonperturbative pion 
distribution amplitude, (j)^{x,Q'^) [2]. The latter can be 
interpreted as the amplitude for the transition of the pion 
with momentum P into two quarks with momenta Px 
and P(l — x). In lowest order pQCD the transition form 
factor is obtained from: 



^ r-0.(x,Q^)+o(a.)+o(%i), (2) 

3 Jo X 



where = 0.131 GeV is the pion decay constant. The 
pion distribution amplitude (DA) plays an important role 
in theoretical descriptions of many hard-scattering QCD 
processes. Since the evolution of (j)T:{x,Q'^) with is 
predicted by pQCD, experimental data on the transition 
form factor can be used to determine its unknown depen- 
dence on X [3-10]. 

The pion transition form factor was measured in the 
CELLO [II] and CLEO [12] experiments in the momen- 
tum transfer ranges 0.7-2.2 GeV^ and 1.6-8.0 GeV^, re- 
spectively. In this paper we study the form factor in the 
range from 4 to 40 GcV^. 



II. THE BABAR DETECTOR AND DATA 
SAMPLES 

We analyze a data sample corresponding to an inte- 
grated luminosity of about 442 fb~^ recorded with the 
BaBAR detector [13] at the PEP-II asymmetric-energy 
storage rings. At PEP-II. 9-GeV electrons collide with 
3.1-GeV positrons to yield a center-of-mass energy of 
10.58 GeV (the T(4S) resonance). Additional data (~ 
43 fb-i) recorded at 10.54 GeV for the purpose of T(4S) 
background studies are included in the present analysis. 

Charged-particlc tracking is provided by a five-layer 
silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH), operating in a 1.5-T axial magnetic field. The 
transverse momentum resolution is 0.47% at 1 GeV/c. 
Energies of photons and electrons are measured with a 
CsI(Tl) electromagnetic calorimeter (EMC) with a res- 
olution of 3% at 1 GeV. Charged-particle identification 
is provided by specific ionization [AE/dx) measurements 
in the SVT and DCH, and by an internally reflecting 
ring-imaging Cherenkov detector (DIRC). 

Signal and background e^e~ — > e"''e~7r°7r'^ processes 
arc simulated with the Monte Carlo (MC) event gen- 
erator GGResRc. It uses the formula for the differen- 
tial cross section from Ref. [1] for 7r° production and the 
BGSM formalism [14] for the two pion final state. Be- 
cause the distribution is peaked near zero, the MC 
events are generated with a restriction on the momentum 
transfer to one of the electrons: ^ —q\ > i GeV^. 
This restriction corresponds to the limit of detector ac- 
ceptance for the tagged electron. The second electron is 
required to have momentum transfer — g| < 0.6 GeV^. 
The experimental criterion providing this restriction for 
data events is described in Sec. III. The pseudoscalar 
form factor is fixed to F{Q) in MC simulation. 

The GGResRc event generator includes next-to- 
leading-order radiative corrections to the Born cross sec- 
tion calculated according to Ref. [15]. In particular, it 
generates extra soft photons emitted by the initial and 
final state electrons. The formulae from Ref. [15] were 
modified to account for the hadron contribution to the 
vacuum polarization diagrams. The maximum energy 
of the extra photon emitted from the initial state is re- 
stricted by the requirement^ E* < 0.05^^, where y^s 
is the e'^e~ center-of-mass (cm.) energy. The gener- 
ated events are subjected to detailed detector simulation 
based on GEANT4 [16], and are reconstructed with the 
software chain used for the experimental data. Varia- 
tions in the detector and beam background conditions 
are taken into account. 

Background events from e+e^ — > qq, where q repre- 
sents a w, d, s or c quark, e^e^ t^t^, and e^e^ — > 
BB are simulated with the JETSET [17], KK2F [18], and 



^ Throughout this paper the asterisk denotes quantities in the 
e+e~ cm. frame. 
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EvtGcn [19] event generators, respectively. 



III. EVENT SELECTION 

At the trigger level candidate events for the process 
under study are selected by the VirtualCompton filter. 
This filter was originally designed to select so-called vir- 
tual Compton scattering (VCS) events used for detector 
calibration. This process corresponds to e^e~ — > e^e~^ 
with the kinematic requirement that one of the final 
state electrons goes along the collision axis, while the 
other electron and the photon are scattered at large an- 
gles. The filter requires that a candidate event contain a 
track with p* / y/s > 0.1 and a cluster in the EMC with 
E* /y/s > 0.1 which is approximately opposite in azimuth 
{\6(j)* — 7r| < 0.1 rad) to this track. Cluster and track 
polar angle acolinearity in the cm. frame is required 
to be greater than 0.1 rad. Finally, the measured miss- 
ing energy in the cm. frame, which should correspond 
to the undetected electron, is compared to a prediction 
based entirely on the directions of the detected particles, 
and the assumption that the missing momentum is di- 
rected along the collision axis: |£^mcas~^prodl/v^ < 0.05. 
For a significant fraction of the e"'"e~ — *■ e~'"e~7r° events 
the trigger cluster algorithm cannot separate the pho- 
tons from TT^ decay, and hence identifies them as a single 
photon. Therefore the VirtualCompton filter has rela- 
tively large efficiency (about 50-80% depending on the 
TT^ energy) for signal events. 

In each event selected by the VirtualCompton fil- 
ter we search for an electron and a 7r° candidate. A 
charged track identified as an electron must originate 
from the interaction point and be in the polar angle range 
0.376 < 6e < 2.450 rad. in the laboratory frame. The 
latter requirement is needed to provide high efficiency 
for the trigger track-finding algorithm and for good elec- 
tron identification. To recover electron energy loss due 
to bremsstrahlung, both internal and in the detector ma- 
terial before the DCH, we look for EMC showers close 
to the electron direction and combine their energies with 
the measured energy of the electron track. The result- 
ing laboratory energy of the electron candidate must be 
greater than 2 GeV. Two photon candidates with en- 
ergies greater than 50 MeV are combined to form a 7r° 
candidate by requiring that their invariant mass be in the 
range 0.06-0.21 GeV/c^ and that their laboratory energy 
sum be greater than 1.5 GeV. 

The main background process, VCS, has a cross section 
several thousand times greater than that for the process 
under study. The VCS photon together with a soft pho- 
ton, for example from beam background, may give an 
invariant mass value close to the tt'^ mass. Such back- 
ground events are effectively rejected by requirements on 
the photon helicity angle (| cos 0^1 < 0.8) and on the 7r° 
cm. polar angle (|cos0*| < 0.8). The photon helicity 
angle 9h is defined as the angle between the decay pho- 
ton momentum in the tt'^ rest frame and the tt*^ direction 
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FIG. 2; The distribution of the cosine of the polar angle of the 
eTr" system momentum in the cm. frame for data (solid his- 
togram) and simulated signal (dotted histogram). Events for 
which I cos 6'*^! > 0.99 (indicated by the arrow) are retained. 



in the laboratory frame. 

The next step is to remove improperly reconstructed 
QED events. We remove events which involve noisy EMC 
channels, events with extra tracks close to the tt" candi- 
date direction, and events with |A0-y^| < 0.025 rad, where 
A9^j is the difference between the laboratory polar an- 
gles of the photons from the 7r° decay. The latter condi- 
tion removes VCS events where the photon converted to 
an e^e^ pair within the DCH volume. It also removes 
about 20% of the signal events, but significantly improves 
(by a factor of about 15) the signal-to-background ratio. 

Two additional event kinematics requirements provide 
further background suppression and improved data to 
MC-simulation correspondence. Figure 2 shows the data 
and MC simulation distributions of the cosine of the po- 
lar angle of the momentum vector of the e7r° system in 
the cm. frame. We require | cos 6**^1 > 0.99. This effec- 
tively limits the value of the momentum transfer to the 
untagged electron ((/I) and guarantees compliance with 
the condition — < 0.6 GeV^ used in MC simulation. 

The emission of extra photons by the electrons involved 
leads to a difference between the measured and actual 
values of Q^. In the case of initial state radiation (ISR) 



)2 

i mcas 



y), where = lEt.jsfs. To restrict 



the energy of the ISR photon we use the parameter 



El 



(3) 



where E*^ and p*^ are the cm. energy and the magni- 
tude of the momentum of the detected eix^ system. In 
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FIG. 3: The r distribution for data events. The shaded his- 
togram shows the background contribution estimated from 
the fit to the two-photon mass distribution (Sec. IV). The ver- 
tical lines indicate the region used to select candidate events 
(-0.025 <r < 0.050). 



FIG. 4: The detection efficiency as a function of the mo- 
mentum transfer squared for events with a tagged electron 
(squares), a tagged positron (triangles), and their sum (cir- 
cles) . 



the ISR case this parameter coincides with defined 
above. The condition r < 0.075 ensures compliance with 
the restriction < 0.1 used in MC simulation. The r 
distribution for data is shown in Fig. 3, where the shaded 
histogram shows the background estimated from the fit 
to the two-photon mass distribution (Sec. IV). We select 
events with —0.025 < r < 0.050 for further analysis. 

The background from e+e^ annihilation into hadrons 
is strongly suppressed by the requirements of electron 
identification, on cos0*^, and on r. An additional two- 
fold suppression of this background is provided by the 
condition that the z-component of the cm. momentum 
of the ctt" system is negative (positive) for events with a 
tagged positron (electron). 

The dependence of the detection efficiency obtained 
from MC simulation is shown in Fig. 4. The detector ac- 
ceptance limits the detection efficiency at small Q^. To 
avoid possible systematics due to data-simulation differ- 
ences near detector edges, we measure the cross section 
and form factor in the region > 4 GeV^. The asym- 
metry of the e+e^ collisions at PEP-II leads to different 
efficiencies for events with electron and positron tags. 
The range from 4 to 7 GeV^ is measured only with 
the positron tag. The decrease of the detection efficiency 
in the region > 10 GeV^ is explained by the decrease 
of the 7r° reconstruction efficiency due to growth of the 
average 7r° energy with Q^. 



IV. FITTING THE TWO-PHOTON MASS 
SPECTRUM 

The two-photon mass spectrum for selected data 
events with 4 < < 40 GeV^ is shown in Fig. 5; for 
> 40 GeV^ we do not see evidence of a 7r° signal over 
background. To determine the number of events contain- 
ing a tt" we perform a binned likelihood fit to the spec- 
trum with a sum of signal and background distributions. 
We describe the signal line shape by a sum of two Fbi 
functions with the same position of their maxima [20]. 
The function Fbi is the convolution of a Gaussian and 
an exponential distribution: 



FBi(a;;Xg,crg, A) = — exp ( ^ — + ^ 



2|A| 



A 



1 — erf 



V2aJA| 



(4) 



The parameters of the 7r° resolution function arc fixed 
from the fit of the mass spectrum obtained for simulated 
signal events weighted to yield the Q'^ dependence ob- 
served in data. The background distribution is described 
either by a linear function in the mass range 0.085-0.185 
GeV/c^ or a second order polynomial in the mass range 
0.06-0.21 GeV/c^. The data mass spectrum is fitted with 
5 (6 for second order polynomial) free parameters: the 
number of signal events, the peak position, the sigma of 
one (narrow) of the i^ei functions (cri), and 2 (3) pa- 
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FIG. 5: The two-photon invariant mass spectrum for data 
events with 4 < < 40 GeV^. The soUd (dotted) curve 
corresponds to the fit with a Unear (quadratic) background 
shape. The dashed curve represents the fitted quadratic back- 
ground. 



rameters for the background. The results of the fits are 
shown in Fig. 5. 

The total number of signal events is about 14000. The 
difference in signal yield between the two background hy- 
potheses is 170 events, while the statistical error on the 
signal yield is 140 events. The difference between the 
peak positions in data and MC simulation is consistent 
with zero. The value of cri is 7.5 MeV/c^ in data and 7.7 
MeV/c^ in simulation, which corresponds to a difference 
of about two standard deviations. 

A similar fitting procedure is applied in each of the 
seventeen intervals indicated in Table I. The param- 
eters of the 7r° resolution function are taken from the fit 
of the mass spectrum for simulated events in the cor- 
responding interval. For the fits to the data, the 
value of the parameter cti is modified to take into ac- 
count the obs erved data-simula tion difference in resolu- 
tion: (Ti \/ — (1-9 MeV)2. The free parameters in 
the data fits are the number of signal events and two 
or three parameters, depending upon the description of 
the background shape. The numbers of signal events ob- 
tained from the fits using a linear background arc listed 
in Table I. The difference between the fits for the two 
background hypotheses is used as an estimate of the sys- 
tematic uncertainty associated with the unknown back- 
ground shape. The two-photon mass spectra and fitted 
curves for three representative intervals are shown in 
Fig. 6. 



V. PEAKING BACKGROUND ESTIMATION 
AND SUBTRACTION 

Background containing true tt'^'s might arise from pro- 
cesses such as beam-gas interaction, e~ annihilation, 
and two-photon processes yielding higher multiplicity fi- 
nal states. 

For beam-gas interactions, the total energy of the de- 
tected electron and 7r° should be less than the beam en- 
ergy. In the energy spectrum of the evr'^-system we do 
not see events with energy less than the beam energy. 
Therefore we conclude that beam-gas background does 
not survive the selection criteria. 

For events due to the signal process with tagged 
positron (electron), the momentum of the detected ctt*^ 
system in the e~ cm. frame has negative (positive) 
z-component, while events resulting from e^e~ annihi- 
lation should be produced symmetrically. Events with 
the wrong sign of the ctt" momentum z-component can 
therefore be used to estimate the background contribu- 
tion from e~ annihilation. The two-photon mass spec- 
trum for such background events is shown in Fig. 7. The 
total number of wrong-sign events is about 3% of the 
selected signal event candidates. The spectrum is fitted 
using a sum of signal and background distributions as 
described in Sec. IV. The fit yields 6 ± 16 7r° events. 
Assuming that the numbers of background events from 

e~ annihilation in the wrong and right-sign data sam- 
ples are approximately the same, we conclude that this 
background does not exceed 0.2% of signal events, and 
so is negligible. Nevertheless we have analyzed simu- 
lated events for the processes e+e" — > BB, e+e" — > qq, 
and e~ — *■ t'^t^ . The number of simulated events for 
each reaction is close to the number of such events pro- 
duced in the experiment. Seven events with right-sign 
z-componcnt of the ctt*' momentum (4 from t'^t~ and 
3 from qq) satisfy the analysis selection criteria, while 
the number of accepted wrong-sign events is four. This 
supports the conclusion that the e^e^ annihilation back- 
ground is negligible. 

The major source of peaking background, of order 10%, 
is two-photon production of two tt^'s. This background 
is clearly seen as a tt'^ peak in the two-photon invariant 
mass spectrum for data events with two extra photons. 
The following procedure is used to estimate the 27r° back- 
ground. We select a clean sample of 27r° events with the 
special selection criteria (described below) and measure 
the distribution for these events {N^-K^i)- Then we 
tune the MC simulation of the e^e~ e"*~e~7r°7r° re- 
action to reproduce the 27r° mass and tt^ angular dis- 
tributions observed in data. Using the MC simulation 
we calculate the ratio (k^) of the numbers of 27r° events 
selected with the standard and special criteria and esti- 
mate the number of 27r° events for each interval that 
satisfy the standard selection criteria as KiN2Tr,i- 

To select 27r'' events we remove the criteria on r and 
cos 91^ and search for events with an extra tt^ . Com- 
binatorial background due to soft false photons is re- 
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FIG. 6: The two-photon invariant mass spectra for data events from three representative intervals. The sohd (dotted) 
curve corresponds to the fit with a Unear (quadratic) background shape. The dashed curve represents the fitted quadratic 
background. 



TABLE I: For each interval, the number of events with tt" obtained from the fit (N^), number of 

background events (A'bkg), total efficiency correction (5totai), number of signal events corrected for data/MC difference and 
resolution effects {Ncoi), and detection efficiency obtained from simulation (e). The quoted errors on A^^ and Ncoi are sta- 
tistical and systematic. For Ncot we quote only Q^-dependent systematic errors. The Q^-independent systematic error is 2.5%. 



interval (GeV^) 




Abkg 




Ncov 


e{%) 


4.0 


-4.5 


1645 ±45 ±4 


176 ±41 


-4.9±1.2 


1503 ±52 ±52 


5.2 


4.5 


- 5.0 


1920±49± 11 


254 ± 54 


-5.5±1.1 


1740 ± 58 ± 70 


9.0 


5.0 


- 5.5 


1646 ±46 ±5 


206 ± 34 


-5.0±1.1 


1551 ±56 ±46 


10.3 


5.5 


-6.0 


1252 ±41 ± 5 


175 ± 30 


-5.5±1.0 


1139 ±50 ±40 


10.7 


6.0 


- 7.0 


1891 ±50 ±2 


271 ±36 


-7.0±1.1 


1760 ± 59 ± 47 


11.5 


7.0 


-8.0 


1229±41 ± 19 


150 ±29 


-7.5±1.0 


1160 ±50 ±44 


13.2 


8.0 


-9.0 


985 ± 38 ± 27 


125 ± 24 


-7.3 ±0.9 


915 ±46 ±46 


15.3 


9.0 


- 10.0 


829 ± 34 ± 8 


59 ±14 


-7.7±1.0 


849 ± 43 ± 23 


17.3 


10.0 


- 11.0 


625 ± 30 ± 18 


47 ±13 


-8.3±1.1 


634 ± 40 ± 30 


18.6 


11.0 


- 12.0 


448 ± 26 ± 3 


27±11 


-8.4±1.0 


484 ± 35 ± 16 


18.7 


12.0 


- 13.5 


405 ± 26 ± 22 


51 ±12 


-8.1±0.9 


381 ±33 ±32 


17.2 


13.5 


- 15.0 


289 ± 22 ± 14 


13±6 


-7.3±1.0 


304 ± 28 ± 20 


16.3 


15.0 


- 17.0 


260 ± 22 ± 5 


14 ±6 


-6.7±1.0 


270 ± 27 ± 11 


15.4 


17.0 


- 20.0 


235 ±21 ± 2 


20±6 


-6.6±1.1 


234 ± 25 ± 10 


13.9 


20.0 


- 25.0 


171 ± 19 ± 11 


5±4 


-6.6 ±1.3 


185 ± 22 ± 14 


11.4 


25.0 


-30.0 


36 ±12 ±2 


1±1 


-6.9±1.5 


36 ± 14 ± 3 


9.2 


30.0 


-40.0 


49 ± 12 ± 2 


2±6 


-6.3±1.8 


53 ±13 ±8 


7.3 



duced by requiring that, in the laboratory, the energy 
of the extra 7r° be greater than 0.2 GeV, and that the 
energies of the decay photons be greater than 50 MeV. 
The mass of the first 7r° must be in the range 0.10-0.17 
GeV/c^. We calculate the parameters cos9* and r for 
the found en^n^ system, and require |cos0*^^| > 0.99 
and —0.025 < r < 0.05. The two-photon invariant mass 
spectrum for the extra tt*^ candidates is shown in Fig. 8. 
The mass spectrum is fitted using a sum of signal and 
background distributions. The signal distribution is ob- 
tained from MC simulation for the e+e^ e+e^Tr^Tr" 
process. The background is described by a second order 
polynomial. This fitting procedure is performed for all 
intervals. To estimate the systematic uncertainty due 



to the unknown background shape we make two fits, one 
with a linear, and one with a quadratic background. The 
difference in tt'' signal size between the fits is taken as a 
measure of systematic uncertainty. 

The 27r° mass spectrum for selected 2tt^ events after 
background subtraction is shown in Fig. 9. The observed 
spectrum differs strongly from the spectrum measured 
in the no-tag mode [21], where the dominant mechanism 
of 27r° production is 77 /2(1270) 7r°7r°. In the 
no-tag mode the /2(1270) meson is produced predomi- 
nantly in the helicity 2 state with angular distribution 
^ sin"' ^TTTT, where the Otttt is the angle between the 
direction and the assumed 77 collision axis in the di- 
pion rest frame. The cos^tttt distribution for selected 2tt° 
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FIG. 7: The two-photon mass spectrum for data with the 
wrong sign of the eTr" cm. momentum component along the 
collision axis. The solid curve is the result of the fit de- 
scribed in the text. The dashed curve represents the fitted 
non-peaking background. 




M^(GeV/c') 



FIG. 8: The invariant mass spectrum of the extra tt candi- 
date in 2-^° data events (points with error bars). The his- 
togram is the result of the fit using a sum of signal and back- 
ground distributions. The dashed curve represents the fitted 
background distribution. 



FIG. 9; The 27r'' invariant mass spectrum for data. 



events after background subtraction is shown in Fig. 10. 
It is seen that our criteria select events with ^tttt near zero 
and strongly suppress /2(1270) production in the helicity- 
2 state. The spectrum in Fig. 9 contains three compo- 
nents: tensor /2(1270), scalar /o(980), and a broad bump 
below 0.8 GeV/c^. We reweight the simulated events to 
reproduce the mass spectrum observed in data. Since the 
mass spectrum may change with Q^, the rcwcighting is 
performed for two intervals (4 < < 10 GeV^ and 
10 < 0^ < 40 GcV^) separately. 

The simulated events are generated with isotropic tt" 
angular distribution in the 27r'^ rest frame. Comparison 
of the simulated cos Ot^tt distribution with the data dis- 
tribution is shown in Fig. 10. We reweight the /2(1270) 
subsample of simulated events so that the total MC sim- 
ulated distribution of Fig. 10 matches the data. Using 
rcwcighted simulated events wc calculate the depen- 
dence of the scale factor Ki which varies from 2.4 at 
Q2 ^ 5 GeV^ to about 1 at > 15 QeV^. The numbers 
of 27r° background events which satisfy our standard se- 
lection criteria are listed in Table I. The fraction of 27r° 
background events in the e7r° data sample changes from 
about 13% for < 10 GeV^ to 6-7% for > 10 GeV^. 

A similar technique is used to search for background 
from the processes e'^e^ —^ e~^e~TT^rj, 77 — > 77 and 
e+e^ —f e'^e^u), uj 77^7. We do not see a clear 77 sig- 
nal in the two-photon mass spectrum, nor do we see an uj 
signal in the 7r''7 mass spectrum; we estimate that these 
backgrounds do not exceed 5% of the 27r° background 
and thus arc negligible. 
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FIG. 10: The distribution of cos 6 
with error bars) and simulated e'^e 
isotropic vr" angular distribution (histogram) 



VI. EFFICIENCY CORRECTION 

The values of the 7*7 tt*^ transition form factor in 
bins of are determined from the ratio of the distri- 
butions from data and MC simulation. The data distri- 
bution must be corrected to account for data-simulation 
difference in detector response: 



for data events (points FIG. 11: The correction to the MC-estimated tt" efficiency as 



N, 



(5) 



where i denotes the interval of under consideration, 
and the Sys are the corrections for the effects discussed 
in Sees. VI A- VXD. 



a function of tt energy in the laboratory frame. 



reconstructed tt^ are selected with our standard criteria 
for the photons and the 7r°, as described in Sec. III. 

The ratio of the reconstruction efficiencies obtained in 
data and in simulation provides a tt*^ efficiency correc- 
tion. This correction. 61 = edata/cMC ~ I7 is shown as a 
fimction of laboratory energy in Fig. 11. The energy 
dependence is well described by a linear function. We 
estimate that the systematic uncertainty associated with 
this correction does not exceed 1%. 

To obtain the correction to the MC-estimated 7r° ef- 
ficiency as a function of we convolve the correction 
energy dependence of Fig. 11 with the 7r° energy spec- 
trum in each interval. The dependence obtained 
is shown in Fig. 12. Only statistical errors are shown; 
the systematic error is estimated to be 1%. 



A. TT reconstruction efficiency 

A possible source of data-simulation difference is 7r° 
loss due to the merging of electromagnetic showers pro- 
duced by the two photons from the 7r° decay, the loss 
of at least one of the decay photons, or the rejection of 
the 7r° because of the selection criteria. The 7r° efficiency 
is studied by using events produced in the ISR process 
e+e" ujj, where uj 7r"''7r~7r° [22]. These events can 
be selected and reconstructed using the measured param- 
eters for only the two charged tracks and the ISR photon. 
Taking the ratio of the number of events with found tt^ to 
the total number of selected e~^e~ —> uj^ events we mea- 
sure the TT^ reconstruction efficiency. The events with 



B. Electron identification efficiency 

The average electron identification (EID) inefficiency 
in the signal MC simulation is about 1%. To estimate 
the data-simulation difference in EID we use VCS events 
which can be selected with negligible background without 
any EID requirements. The EID efficiency is determined 
as the ratio of the number of events with an identified 
electron to the total number of VCS events. The ratio of 
the efficiencies obtained from data and simulation gives 
the efficiency correction. We determine the correction as 
a function of the electron energy and polar angle and con- 
volve this function with the electron energy and angular 
distributions for the process under study. The resulting 
dependence of the efficiency correction is shown in 
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FIG. 12: The dependence of the correction to the MC- 
estimated efficiency for 7r°. 
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FIG. 13: The correction to the MC-estimated EID efficiency 
as a function of Q^. 



Fig. 13. 



C. Trigger efficiency 

With the available statistics and the trigger configu- 
ration used, we cannot determine the trigger efficiency 



FIG. 14: The dependence of the correction for trigger effi- 
ciency for events with a tagged electron (triangles) or positron 
(squares), and for all events (circles). 



for the process under study by using the data. However 
the trigger efficiency can be measured for the VCS pro- 
cess, which has a much larger cross section. The VCS 
events allow the determination of the part of the trigger 
inefficiency related to the trigger track-finding algorithm. 
The remaining trigger inefficiency, which is related to 
the ability of the trigger cluster algorithm to separate 
nearby photons from tt*' decay, depends strongly on tt^ 
energy. Therefore the data-simulation difference can be 
estimated from comparison of the tt*^ energy spectra in 
data and simulation. 

The VCS events are selected with the criteria described 
in Sec. Ill after applying the 7r° requirements to the VCS 
photon. These events must satisfy a second trigger line 
that selects VCS events with an efficiency close to 100%, 
but that is prescaled by a factor of 1000. The trigger 
efficiency is determined as the fraction of selected events 
which pass the VirtualCompton filter. The ratio of the 
efficiencies obtained from data and simulation gives the 
efficiency correction. We find that the trigger efficiency 
depends strongly on electron scattering angle. When 9e 
changes from 0.376 to 0.317 rad, the efficiency falls from 
70% to 30% and the efficiency correction increases from 
10% to 30%. For this reason the events with Be < 0.376 
rad were removed from the analysis data sample. Since 
the angular and energy distributions for the VCS and 
e'^e~ — > e"^e~7r'^ processes are different, we determine 
the correction as a function of cos 9^ and cos^e, sepa- 
rately for tagged electrons and positrons, and convolve 
this function with the cos and cos 9^ distributions for 
the process under study. The resulting dependence 
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FIG. 15: The trigger inefficiency as a function of 7r° energy 
determined directly from e^e~ — > e'^e~Tv° simulation (circles) 
and from simulated VCS events (squares). 



FIG. 16: The spectrum for simulated signal events. The 
shading represents the contributions from the different vr" en- 
ergy ranges (in GeV) indicated by the key. 



of the efficiency correction is shown in Fig. 14. The cor- 
rections for events with a tagged electron or positron are 
also shown. The correction for tagged positron events is 
about —2% and flat. For events with a tagged electron, 
the graph begins at = 7 GcV^. The electron correc- 
tion changes from —8% at = 7 GeV^ to about —1.5% 
at Q2 = 20 GcV^ and higher. 

The trigger inefficiency determined directly from 
e^e^ — > e^e~7r*^ simulation is compared to that calcu- 
lated using simulated VCS events in Fig. 15. The dis- 
crepancy between the inefficiencies is 3-4% for ir^ ener- 
gies higher than 3 GeV, but increases to 30% for i?^ < 2 
GeV. Figure 16 shows the spectra for simulated sig- 
nal pions of different energies. For each interval the 
energy spectra for data and simulation should be identi- 
cal. Any difference provides a measure of the quality of 
the trigger efficiency simulation. Using the fitting proce- 
dure described in Sec. IV, we determine the numbers of 
signal events in ten intervals for six tt'^ energy ranges 
(53 measurements, excluding cells with no events). We 
subtract the 27r° background and normalize the energy 
spectrum in each interval so that its integral is unity. 
The same procedure is applied to simulated spectra af- 
ter introducing the efficiency corrections for 7r° loss, EID 
and trigger inefficiency. The comparison of the normal- 
ized data and simulated spectra gives x^/ndf — 42.4/43 
{ndf = number of degrees of freedom). The 7r° energy 
spectrum summed over the intervals from 4 to 12 
GcV^ is shown in Fig. 17. The simulated spectrum is 
the sum of the spectra normalized to the number of data 
events in each interval. The shaded boxes represent 



the uncertainties in the efficiency corrections. The ratio 
of the data and simulated spectra agrees with unity with 

/ndf =■ 5.9/5. Since the spectra for data and sim- 
ulation are in reasonable agreement, we conclude that 
the simulation reproduces the discrepancy in trigger in- 
efficiency, and so there is no need to introduce an extra 
efficiency correction. However, we introduce an extra sys- 
tematic uncertainty due to trigger inefficiency, which is 
conservatively estimated to be 2%, i.e., half of the dif- 
ference between the VCS and e^e~ e^e~7r'' trigger 
inefficiencies for high energy tt^'s (see Fig. 15). 

Since the energy and angular distributions, and the 
trigger efficiency correction are very different for events 
with a tagged electron or a tagged positron, it is in- 
teresting to compare the e'^e~ e^e~7r*' differential 
cross sections measured for electron-tagged and positron- 
tagged events. To do this we subtract the 27r° back- 
ground and apply the efficiency correction separately to 
events with tagged electron and positron. The ratio of 
the cross sections is then calculated as the double ratio 

(A'"data/^Mc)electrons/(A^data/-/VMc)positrons (Fig. 18) and 

is found to be in reasonable agreement with unity. 



D. Requirements on r and cos^e^ 

To estimate possible systematic uncertainty due to 
the requirement —0.025 < r < 0.05 we study events 
in the range 0.05 < ?' < 0.075 (see Eq. (3) and 
Fig. 3). Assuming that the efficiency corrections are 
the same for events from the two r ranges, and sub- 
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FIG. 17: The vr" energy spectrum for data (points with error 
bars) and signal simulation (histogram). The shaded boxes 
represent the uncertainties associated with the efficiency cor- 
rection procedure. 
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FIG. 19: The dependence of the total efficiency correction. 



trading 27r*' background, we determine the double ratio 

(-^data/-^Mc)o.05<r<0.075/(-^data/iVMc)-0.02.5<r<0.05 aS a 

function of Q^. The ratio is consistent with unity 
{x^/ndf = 9.3/15), so we conclude that the MC simu- 
lation reproduces the shape of the r distribution. 

We also study the effect of the | cos0*^| > 0.99 crite- 
rion by changing the value to 0.98 and 0.95. The ratio 
of the numbers of events with 0.98 < |cos0*^| < 0.99 
and I cos 6**^ I > 0.99 is found to be 0.013 ± 0.003 in data 
and 0.0074 ± 0.0004 in simulation. The corresponding 
values for 0.95 < | cos 6**^1 < 0.99 are 0.018 ± 0.002 and 
0.0103 ± 0.0005. Since the observed data-simulation dif- 
ference does not exceed 1%, we do not introduce an ef- 
ficiency correction, but consider this difference (1%) as 
a measure of systematic uncertainty due to the cos0*^ 
criterion. 

The total efficiency correction as a function of is 
shown in Fig. 19, where the error bars are of statistical 
origin. The systematic uncertainty, which is independent 
of Q^, is 2.5% and takes into account the uncertainties in 
the determination of tt" loss (1%) and trigger inefficiency 
(2%), and the uncertainty due to the jcos^*^! > 0.99 
requirement (1%). The values of the total efficiency cor- 
rection and their statistical errors are listed in Table I. 



FIG. 18: The dependence of the double ratio defined 
in the text to compare the cross section for selected signal 
events with a tagged electron to that for events with a tagged 
positron. 



VII. CROSS SECTION AND FORM FACTOR 



The Born differential cross section for e^e 



e ' e n" 



is calculated as 



da ^ iVcor/AQ' 
dQ2 eRL 



(6) 
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20: The dependence of the radiative correction fac- 



where iVcor is the number of signal events corrected 
for data-simulation difference and resolution effects (Ta- 
ble I), AQ^ is the relevant interval, L is the total 
integrated luminosity (442 fb~^), e is the detection effi- 
ciency as a function of , and i? is a radiative correction 
factor accounting for distortion of the spectrum due 
to the emission of photons from the initial state particles 
and for vacuum polarization effects. The detection effi- 
ciency is obtained from simulation. Its dependence is 
shown in Fig. 4 and listed in Table I. 

The radiative correction factor is determined using 
generator-only simulation. The spectrum is generated 
using only the Born amplitude for the e^e" — *■ e+e^Tr" 
process, and then again using a model with radiative cor- 
rections included. The dependence of the radiative 
correction factor, evaluated as the ratio of the second 
spectrum to the first, is shown in Fig. 20. The depen- 
dence is fitted by the function a/ {1 + bQ'^). The accuracy 
of the radiative correction calculation is estimated to be 
1% [15]. Note that the value of R depends on the require- 
ment on the extra photon energy. The dependence 
obtained corresponds to the criterion r = 2E*/^ < 0.1 
imposed in the simulation. 

The corrected mass spectrum (A^cor) is obtained from 
the measured spectrum (A^rec) by dividing by the effi- 
ciency correction factor [see Eq.(5)] and unfolding the ef- 
fect of resolution. Using MC simulation, a migration 
matrix A is obtained, which represents the probability 
that an event with true in interval j is reconstructed 
in interval i: 



TABLE II: The 

for the interval (Q^), the e^e 
{da/dQ^(Q'^)), and the product of the 77* — > tt" transition 
form factor F{Q^) and Q^. The quoted errors are statistical 
and systematic for the cross section, and combined for the 
form factor. In the table we quote only Q^-dependent sys- 
tematic errors. The Q^-independent systematic error is 3% 
for the cross section and 2.3% for the form factor. 
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(7) 



In the case of extra photon emission, Q^i-^g is calculated 
as —{p — p' — k)^, where k is the photon four-momentum; 
e and R in Eq.(6) arc functions of Qtruc- The reso- 
lution varies from about 0.05 GeV^ at = 5 GeV^ to 
0.25 GeV^ at = 25 GeV^. As the chosen inter- 
val size significantly exceeds the resolution for all , the 
migration matrix is nearly diagonal, with diagonal values 
0.9, and ncxt-to-diagonal values ^ 0.05. We unfold the 
spectrum by applying the inverse of the migration ma- 
trix to the measured spectrum. The procedure changes 
the shape of the distribution insignificantly, but in- 
creases the errors (by about 20%) and their correlations. 
The corrected spectrum is listed in Tabic I. 

The values of the differential cross section are listed 
in Table II. The quoted errors are statistical and sys- 
tematic. The latter includes only Q^-dependent errors: 
the systematic uncertainty in the number of signal events 
and the statistical error in the detection efficiency deter- 
mined from MC simulation. The Q^-independent sys- 
tematic error is equal to 3% and includes the systematic 
uncertainties in the efficiency correction (2.5%) and in 
the radiative correction factor (1%), and the uncertainty 
in the integrated luminosity (1%). The measured differ- 
ential cross section at the Born level is shown in Fig. 21, 
together with CLEO data [12] for > 4 GeV^ 

Because of the strong nonlinear dependence of the 
cross section on , the effective value of correspond- 
ing to the measured cross section differs from the center of 
the interval. We parametrize the measured cross sec- 
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FIG. 21: The e"*"e~ e^e~-K^ differential cross section ob- 
tained in this experiment compared to that from the CLEO 
experiment [12]. 



tion by a smooth function, reweight the distribution 
in simulation to be consistent with data, and calculate 
the weighted average (Q^) for each mass interval. The 
values of are listed in Table II. 

Since the requirement on cos 9*^ limits the momentum 
transfer to the untagged electron, we measure the cross 
section for the restricted (j| range IqII < imax- The value 
of Omax is determined from the q| dependence of the de- 
tection efficiency (e(aniax) = 50%) and is equal to 0.18 
GeV2. 

To extract the transition form factor we compare the 
measured and the calculated values of the cross section. 
The simulation uses a constant form factor Fmc- There- 
fore the measured form factor is determined from 



(dcr/dQ2)data 



(da/dQ2)MC 



p2 

2V ^MC- 



(8) 



The calculated cross section (dCT/dQ^)MC has a model- 
dependent uncertainty due to the unknown dependence 
on the momentum transfer to the untagged electron. 
We use a gl-independent form factor, which corresponds 
to the QCD-inspired model F{q1,ql) cx l/{ql + q^) ~ 
l/qf [23]. Using the vector dominance model with the 
form factor -F'(ql) cx; 1/(1 — q'^/'n^'j,), where rrip is p me- 
son mass, leads to a decrease of the cross section by 
3.5%. This difference is considered to be an estimate of 
model uncertainty due to the unknown dependence. 
However, it should be noted that this estimate depends 
strongly on the limit on gl. The value of 3.5% is obtained 
with I^jI < 0.18 GeV^. For a less stringent constraint, 
for example |g|| < 0.6 GeV^, the difference between the 



FIG. 22: The 77* tt'^ transition form factor multiplied by 
Q^. The dashed line indicates the asymptotic limit for the 
form factor. The dotted curve shows the interpolation given 
by Eq.(9). 



calculated cross sections reaches 7.5%. 

The values of the form factor obtained, represented in 
the form Q'^\F{Q'^)\, are listed in Table II and shown 
in Fig. 22. For the form factor we quote the combined 
error, for which the statistical and Q^-dcpendcnt sys- 
tematic uncertainties are added in quadrature. The Q^- 
independent systematic error is 2.3%, and includes the 
uncertainty on the measured differential cross section, 
and the model-dependent uncertainty due to the un- 
known g| dependence. 



VIII. CONCLUSIONS 



We have studied the e+e 



e~^e tt" reaction in the 



single tag mode and measured the differential cross sec- 
tion (da/dQ^) and the 77* 7r° transition form factor 
F{Q^) for the momentum transfer range from 4 to 40 
GeV^. For the latter, the comparison of our results with 
previous measurements [11, 12] is shown in Fig. 22. In 
the range from 4 to 9 GeV^ our results are in rea- 
sonable agreement with the measurements by the CLEO 
collaboration [12], but have significantly better precision. 
We also significantly extend the region over which the 
form factor is measured. 

To effectively describe the dependence of the form 
factor in the range 4-40 GeV^, we fit the function 



Q^\Fm\^A 



10 GeV^ 



/3 



(9) 
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S. V. Mikhailov and N. G. Stefanis using the light-cone 
sum rule method [4, 25] at next-to- leading order (NLO) 
pQCD; the power correction due to the twist-4 contribu- 
tion [25] was also taken into account. Their results are 
shown for the Chernyak-Zhitnitsky DA (CZ) [26], the 
asymptotic DA (ASY) [27], and the DA derived from 
QCD sum rules with non-local condensates (BMS) [28]. 
For all three DAs the dependence is almost flat for 
> 10 GeV^, whereas the data show significant growth 
of the form factor between 8 and 20 GeV^ . This indicates 
that the NLO pQCD approximation with twist-4 power 
correction, which has been widely used for the description 
of the form-factor measurements by the CLEO collabo- 
ration [12], is inadequate for less than ~ 15 GeV^. In 
the range from 20 to 40 GcV^, where uncertainties 
due to higher order pQCD and power corrections are ex- 
pected to be smaller, our data lie above the asymptotic 
limit and arc consistent with the CZ model. 
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FIG. 23: The 77* tt transition form factor multiplied 
by Q^. The dashed line indicates the asymptotic limit for 
the form factor. The solid and dotted lines show the predic- 
tions for the form factor [8] for the CZ [26] and asymptotic 
(ASY) [27] models for the pion distribution amplitude, re- 
spectively. The shaded band represents the prediction for the 
BMS [28] pion DA model. 



to our data. The values obtained for the parameters are 
A = 0.182 ± 0.002 GeV, and /3 = 0.25 ± 0.02. The fit 
result is shown in Fig. 22 by the dotted curve. The ef- 
fective dependence of the form factor (~ dif- 
fers significantly from the leading order pQCD prediction 
(~ 1/Q^) (see Eq.(2)), demonstrating the importance of 
higher-order pQCD and power corrections in the re- 
gion under study. 

The horizontal dashed line in Fig. 22 indicates the 
asymptotic limit Q'^F{Q'^) = V2f^ « 0.185 GeV for 
^ 00, predicted by pQCD [2]. The measured form 
factor exceeds the limit for > 10 GeV'^. This con- 
tradicts most models for the pion distribution amplitude 
(see, e.g., Ref. [24] and references therein), which give 
form factors approaching the asymptotic limit from be- 
low. 

The comparison of the form factor data to the pre- 
dictions of some theoretical models is shown in Fig. 23. 
The calculation of [8] was performed by A. P. Bakulev, 
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